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I. INTRODUCTION 
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In recent years, there is growing interest in higher-dimensional cosmological models inspired by the recent devel- 
. opments in string theory. One sort of such models are branes which are submanifolds embedded into the higher- 
^ dimensional spacetime. They play an important (even essential) role, for instance, in a description of the confinement 
^ of the nongravitational interactions and/or stabilization of the extra dimensions. Braneworld models, in which our 
' universe is exactly one of four-dimensional branes, are the central issues in this research field. Much efforts to reveal 
~l ] cosmology on the brane have been done in the context of five-dimensional spacetime, especially after the stimulating 
• ' proposals by Randall and Sundrum [ij (see, e.g., But of course there is no particular reason to restrict the 

number of dimensions to be five. In recent years, the focus has turned to higher-dimensional braneworld models. 

It is still unclear how the gravity and cosmology on the brane embedded in the higher-dimensional spacetime 
behave. In higher dimensions, there would be many kinds of classes of possible braneworld models. Here we focus on 
the generic feature of higher-dimensional braneworld models and from this we deduce the subject to be revealed in 
this paper. 

It is commonly assumed that a brane is an infinitely thin object. In five dimensions, this is a good approximation 
' as long as one is interested in the scales larger than the brane thickness. Then, the cosmology on the brane in five 
' dimensions, i.e., dynamics of the brane in one bulk spacetime, is uniquely determined by the Israel junction conditions 
once the matter on the brane is specified. But, in higher dimensions the situation is quite different. It happens because 
self-gravity of an infinitesimally thin three-brane in a higher-dimensional spacetime develops a severe singularity and 
one cannot put any kind of matter on the brane if the number of codimensions is larger than 2. Also in the case of 
a codimension-two brane, one cannot put any kind of matter other than the pure tension. The cosmology on such 
a brane must be investigated under some prescription of the singular structure of the brane. The well-motivated 
prescription is to regularize the brane by taking its microscopic structure into consideration. However, in the case of 
six-dimensional models, the cosmology on the brane strongly depends on the way of regularization. This implies that 
there is no unique brane cosmology in the case of higher codimensions. 

Thus, in this paper we take a different point of view. We shall take the brane thickness into account from the 
beginning and look for the regular braneworld solutions with finite thickness, called thick brane solutions, in a field 
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theory model. The thickness of the brane may be very close to the length scale of the quantum gravitational theory, 
e.g., the string length scale. As we mentioned, the brane thickness must be more essential in higher-dimensional 
spacetimes than in the five-dimensional one. Thus, we look for the thick brane solutions in higher dimensions as well 
as in five dimensions. Bearing the cosmological applications, we focus on the thick brane solutions which have de 
Sitter geometry in the ordinary four-dimensions. 

The properties of thick branes in five-dimensions, especially localization of gravity, have been discussed in Ref. 
In five dimensions, the thick brane solutions in five-dimensional spacetime have been explored in various literatures 
(see, e.g., fi|). In this paper, we consider thick brane solutions supported by two interacting scalar fields (j) and x, 
whose potential is given by 

x) = - + ^(.X'- mlf \,cj^\^ + A. (1) 

A stronger interaction between two scalar fields develops the local metastable vacua, which may be able to support 
our brane universe, other than the global one. Several works on thick brane solutions in such a setup have already 
been performed. Five-dimensional Minkowski brane solutions were derived in [^0]. Bearing cosmological applications 
in mind, our goal in this work is to find thick de Sitter brane solutions. Thick brane solutions in higher-dimensional 
spacetimes with two interacting scalar field were also investigated in [7| (see, e.g., for other kinds of models). We 
shall give de Sitter versions of these solutions. 

In Ref. Q there are some arguments in favor of the fact that the scalar fields, used in [7], are a quantum nonpertur- 
bative condensate of a SU(3) gauge field. Briefly these arguments consist of the following: components of the SU(3) 
gauge field can be divided in some natural way onto two parts. The first group contains those components which 
belong to a subgroup SU{2) e SU{3). The remaining components belong to a factor space SU{3)/ SU{2). Following 
Heisenberg's idea [10l| about the non-perturbative quantization of a nonlinear spinor field the nonperturbative quan- 
tization for the SU(3) gauge field is being carried out as foUowings. It is supposed that two-point Green functions 
can be expressed via the scalar fields. The first field (p describes two-point Green functions for SU(2) components of a 
gauge potential, and the second scalar field x describes two-point Green functions for SU{3)/ SU{2) components of the 
gauge potential. It is supposed further that four-point Green functions can be obtained as some bilinear combination 
of the two-point Green functions. Consequently, the Lagrangian of the SU(3) gauge field takes the form It allows 
us to regard such sort of thick brane solutions as some defect in a spacetime filled by a condensate of the gauge field 
living in the bulk. 

This paper is organized as follows. In Sec. II, we introduce the basic field theory model to find thick brane solutions. 
In Sec. Ill, we present thick brane solutions in five-, six- and seven-dimensional models. In Sec. IV, we investigate 
perturbations in the five-dimensional model and then give speculations about stability of higher-dimensional models. 
Secion V contains a brief summary. 



II. BASIC THEORY 



Our interest is in the D = [A + ?i)-dimensional Einstein-scalar theory whose action is given by 
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where M„+4 is the gravitational energy scale in the spacetime with n extra dimensions. The potential energy V is 
defined by Eq ([T]). e = -|-l(— 1) corresponds to the case of normal (phantom) scalar fields and A is an arbitrary constant. 
As we mentioned in the introduction, in Ref. [9], it was argued that such a theory composed of two interacting scalar 
fields coupled to the remaining U{1) degrees of freedom appears as a result of the gauge condensation in the original 
SU{3) gauge theory. In this case, the scalar fields and x correspond to the vacuum expectation values of the SU (2) 
and SU{3)/{SU{2) x U{1)) gauge sectors, respectively. 

We assume that the generalized D-dimensional metric has the static form [TTI| 

ds^ = a'^{r)ja^{x'')dx°'dx'^ - \{r){dr^ + r^dni_^) , (3) 

where d^^_i is the solid angle for the (n — 1) sphere, is the metric of the four-dimensional de Sitter space whose 
scalar curvature is given by R^v['y\ = 2tH'^^^u- The corresponding Einstein equations are given by: 
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The scalar field equations are 
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The extreraum of the bulk potential are located at the following four points 
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We call these points 0, 1, 2 and 3, respectively, with correspoding potentials Vi {i — 0, 1, 2, 3). Clearly, max(V2, V3) < 
Vi. We also find 
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The potential ([T]) has two global minimums 3 and two local minima 2 at values of the parameters Ai , A2 used in 
this paper. The conditions for existence of local minima are: Ai > Q,ra\ > A2to|/(2A3), and for global minima: 
A2 > 0, m| > Aim^/ (2A3). From comparison of the values of the potential in global and local minima we have (taking 
into account that Vioc = at Vb = A2m|/4, see below for the six- and seven-dimensional cases) Vgi < Vioc Vgi < 0, 
i.e. A2to| > Ximf. Besides two local and two global minima, four unstable saddle points (0) and the local maximum 
(1) exist. In the next sections, we look for thick brane solutions starting and finishing in one of local minima 2 of the 
potential HI). 

For our purpose, in the following sections the bulk coordinate and scalar field variables are rescaled to make 



them dimensionless as 
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parameters are also rescaled as V 
J = 1,2. 

Then, the boundary conditions at r — + are given by a'(0) — Q, (j)' — Q and x' = 0- And from the constraint 
relation of Einstein equations Eq. ([5]) , we find the boundary value of scale factor 
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At the asymptotic infinity, the scalar fields approach the local minimum where 0'(oo) — 'x' {00) — 0. 
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III. THICK DE SITTER BRANE SOLUTIONS 
A. The five-dimensional model 



By substituting n = 1 and X{r) ~ 1 in the metric ansatz Eq. ([3]), we have 
Then, the Einstein equations (H]), ([5]) and ^ now read 
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Similarly, the scalar field equations read 



— ^ 
a 

4a' 



X 



-X = X 



2A3x' + Ai(02-m2) 
2A3(^2 + A2(x'-m2) 



(14) 



We solve these equations under the boundary conditions shown at the end of the previous section. Wc will look for 
the thick de Sitter brane solution with the asymptotically anti- de Sitter bulk. 



B. The six-dimensional model 

By substituting n = 2 into our metric ansatz given by ((3]), the general form of the six-dimensional metric is given 

ds^ = airf^^.^dx^'dx" - \{r){dr'^ + r^dd^). (15) 

To find thick brane solutions, we solve the coupled Einstein-scalar system given in Sec. III. By substituting n — 2 
into the Einstein and scalar field equations and combining them appropriately, we find the equations 
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We solve these equations under the boundary conditions shown at the end of the previous section. We will look for 
the thick de Sitter brane solution with the asymptotically flat bulk. 



The seven-dimensional model 



Similarly to the case of the six-dimensional model, we consider the seven-dimensional model. The bulk metric is 
given by Eq. ([3]) with n — S. Then, by combining a set of equations appropriately, the equations of motion can be 
obtained as follows 
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FIG. 1: The scalar field configurations <jf> and X are shown FIG. 2: The metric function A(r) is shown as functions 
as functions of the dimensionless r. of dimensionless r. 

The constraint equation is given by 




As in the five and six-dimensional models, we will solve them under the boundary conditions shown at the end of the 
previous section. We will look for the thick de Sitter brane solution with asymptotically flat bulk. 

D. Numerics and solutions 

Let us present now the numerical examples of thick de Sitter brane solutions for the case of phantom scalar field 
e = — 1. We solved the coupled Einstein-scalar system by the iteration method. The detailed description of the 
method can be found, for example, in Ref. 3- The essence of this procedure is the following: on the first step one can 
solve the equation for the scalar field 4> with some arbitrary selected function x looking for a regular solution existing 
only at some value of the parameter mi. On this step the infiuence of gravitation is not taken into account. Then one 
inserts this solution for the function (j) into the equation for x Emd searches for a value of the parameter TO2 yielding a 
regular solution. This procedure is repeated (three times is usually enough) for obtaining acceptable convergence of 
values of the parameters mi, TO2. The obtained functions (j), x ^^re inserting into the Einstein equations. Equation ([5]), 
which is the constraint equation, is using for specifying the boundary conditions (see above). The obtained solutions 
for the metric functions A, a are inserted into the equations for the scalar fields, and they are solved again for a search 
of eigenvalues of the parameters toi, m2 with an account of gravitation. This procedure is repeated as many times as 
it is necessary for obtaining of acceptable convergence of values of the parameters mi, TO2. 

We set the parameters as Ai = 0.1, A2 = 1.0, A3 — 1.0. For the five-dimensional case we have chosen ^(0) = 1.0, 
x(0) — \/0.6 and A = —3.6. For these parameters we solved the nonlinear eigenvalue problem for mi and m2. 
After 11 iterations, the eigenvalues mi k, 1.989512 and m2 ~ 1.964764 were found. In Fig. 1-4, we showed our 
numerical solutions. The spacetime is asymptotically an anti-de Sitter one: the asymptotic value of the potential ([1]) 
eVoo = e{X2m\/A -I- A) < plays the role of a negative cosmological constant. 

For the six-dimensional case we take the boundary condition 4>{fi) — a/3 « 1.73205 and x(0) = a/O^ ~ 0.774597. 
We also choose the bulk cosmological constant as A = — A2m2/4, in order for our solutions to be asymptotically flat 
in the bulk. After ten iterations, we find the eigenvalues mi « 2.3590 and m2 « 3.0599. In Fig. 1-4 we showed our 
numerical solution for a given set of parameters. 

As for the case of the seven-dimensional case we also choose the bulk cosmological constant as A = — A2m|/4. After 
ten iterations, we find the eigenvalues mi k, 2.24633 and m2 ~ 3.11911. The obtained solutions for this case are 
presented in Figs. 1-4. 
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FIG. 3: The function a ja is shown as functions of di- 
mensionless r. 
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FIG. 4: The energy density = e + x") + V\ 

is shown as functions of dimensionless r. 



IV. STABILITY 



For stability analysis, let us write down the scalar field and Einstein equations, by keeping the time-dependence of 
the metric and field functions. Using the same rescaling as above and introducing new metric functions 

one can get the following (r, t) and (r, r) components of the Einstein equations 
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Here the dot refers to derivative with respect to dimensionless time t = Rt. The scalar filed equations are: 
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Now we search for regular perturbed solutions of the above equations in the form 

a = ao{r) + Sa{T,r), jS ^ l3o{r) + SP{T,r), 
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where the subscript refers to static solutions of the equations and 6a{T^ r)^ SP{t, r) , 6(f){T, r) , Sx{t, r) are perturba- 
tions. Then we will have from (flQl) and (l20l) : 
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and from (ED) and (1221) follows 
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One can see that Eqs. ((25)) - ([27)) contain similar terms like (^Sy + 35y^ , where Sy is one of fimctions 

Sa{T,r),SP{T,r),6ip{T,r),Sx{T,r). Such term implies that one can search for a solution in the form: 

Sy ~ yi{r)e^i'^ cosujT, where oj = — | > 0. It implies that o^o > |- Inserting the last expression into ((24)) - 
(|27ll . one can get 
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So we have four usual differential equations for perturbations. For existence of stable solutions we need to provide 
positiveness of an eigenvalue Wq. To do this we will search for finite solutions of above equations. 

For the five-dimensional model, by setting (3i — Q (and /3o = 0) and substituting n — \, the perturbations equations 
are simply reduced to 



<^i" + 4ao'0i' - [2A3X0 + y'/'o" + Ai (3^2 - ml) ] 0i - (y^o'xo' + 4A3(/)oXo)xi + ^oe"'""'/' 



0, 



Xi" + 4ao'xi' - [2A30^ -I- ^xo" + A2 {ixl " ml) ] xi - (y 0o'xo' + 4A30oXo) 0i + ^oe^'^'Xi = 0, (32) 
where we have used Eq. ([28]) 

3a'i = -e((/>o'?!>i +X0X1) 

to eliminate ai' . We solved numerically the Eqs. ([32]) . using the background solutions obtained in Sec. Ill, with the 
boundary conditions 



0i(O) = l.O, 0UO) = O.O, xi(0) = -1.0, x'i(0) = 0.0, 0i(r->oo) -.0, xi('^ ->«))-. 0. 
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In the case of the thick Minkowski brane solutions supported by the phantom scalar fields, it has been shown that 
these solutions are stable @. In our case, the perturbation equations also have regular solutions. Also, we have 
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FIG. 5: The perturbed scalar field configurations <^i (solid 
curve) and xi for '^o ~ 13.179575 (dashed-dotted curve) are 
shown as functions of r. As the background solutions, we used 
the one shown in Figs. 1-3. 




confirmed that there are no solutions with Uq < 0, i.e. unstable solutions are absent. In Fig 5, we showed an example 
of numerical solutions for perturbations with loq w 13.179575. We conclude that our thick de Sitter brane solutions 
supported by the phantom scalar fields are stable, irrespectively of the value of the Hubble parameter H. As a special 
limit if — > 0, we may recover the result for the Minkowski brane solutions. 

It is also reasonable to expect that our phantom thick brane solutions in six- and seven-dimensional spacetime are 
also completely stable, since the spacetime structures of higher-dimensional solutions are essentially very similar to 
that of the five-dimensional solution. 



V. SUMMARY AND DISCUSSIONS 



In this paper, we have presented five-, six-, and seven-dimensional thick de Sitter brane solutions supported by two 
interacting (phantom) scalar fields. The special form of the potential ([1]) allows us to find regular solutions with finite 
energy density. It was shown that asymptotically there exist anti-de Sitter spacetime for the five-dimensional case 
and flat spacetime for the six- and seven-dimensional cases. 

Also we have investigated the stability of our thick brane solutions. In the five-dimensional thick de Sitter brane 
solutions, we have shown the presence of decaying solutions for perturbations. We also could not find any mode with 
negative eigenvalue. This, of course, implies that our solutions in five-dimensional spacetime are stable. It is also 
quite reasonable to expect that our thick brane solutions in higher dimensions are also stable since the spacetime and 
vacuum structures are essentially the same as in the case of five spacetime dimensions. 
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